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Abstract. We study structures of general relativistic compact stars with exotic matter. Our study 
is based on axisymetric and stationary formalism including purely toroidal magnetic field. We also 
study the finite size effects of quark-hadron mixed phase on structures of magnetars. For hybrid 
stars, we find a characteristic distribution of magnetic field, which has a discontinuity originated 
in the quark-hadron mixed phase. These distributions of magnetic field will change astrophysical 
phenomena, such as cooling processes. 
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INTRODUCTION 

It was presented that quark matter may exist in compact stars 010, Si- So far, there has 
been extensive work devoted to studying the effects of quark matter on astrophysical 
phenomena; the gravitational wave radiation [0, S @], cooling processes [0, H H, [Toll , 
neutrino emissions iflll 121. and rotational frequencies lfl3ll . the maximum energy re- 



lease by conversions from neutron stars to quark/hybrid stars 11141 11511. etc.. However, 
uncertainties of equation of states (EOSs) have been still left. 

For such studies on compact stars, general relativistic effects are fundamentally im- 
portant, since baryon density is comparable to pressure, poc 2 ~ P, hence the gravitational 
forth is strong. Moreover, strong magnetic field may change hydrostatic equilibriums of 
compact stars. 

In this paper, we study the structures of general relativistic hybrid stars with purely 
toroidal magnetic field. To calculate magnetized compact stars, we adopt the scheme 
based on axisymetric and stationary formalism including purely toroidal magnetic 
field[Q3]. A star with pure toroidal magnetic field is unstable, however it becomes an- 
other stable star in which toroidal magnetic field is dominant in the dynamical simu- 
lation fl7|l . Moreover, Heger et al. have suggest that the toroidal magnetic field may 
dominate 10 5 times lager than the poloidal magnetic field at the last stage of the main 
sequence |fl8|l . 

The organization of this paper is as follows. Adopted EOSs are briefly discussed in 
Sec. 2 with equilibriums of magnetized rotating compact stars. In Sec. 3, we show our 
nuemirical results. In Sec. 4, we discuss the consequence of our calculations. 



INPUT PHYSICS 



Equation of state 

The hardness of EOSs is an important ingredient for determining the equilibrium 
configurations as mentioned in Sec. 1. Though proto-neutron stars left after supernova 
explosions are very hot, T ~ 50 MeV, they cool down to cold neutron stars (T ~ 
lMeV) in some tens of seconds Therefore we assume that the temperature for 
hydrostatic compact stars is zero, since EOSs at such low temperature (T ~ lMeV) 
show almost same stiffness as at zero-temperature. It means the temperature of compact 
stars is always much smaller than the typical chemical potentials. Moreover, the loss 
of neutrinos means that the chemical potentials of the neutrinos may be set to zero. 
Thus, we impose the barotropic condition of the EOSs (P = P(e)) by assuming zero- 
temperature, zero-neutrino fraction, and beta-equilibrium. 

Our theoretical framework for the hadronic phase of matter is the nonrelativistic 
Brueckner-Hartree-Fock approach including hyperons such as E~ and A |[20ll . 

For quark phase, we adopt the MIT bag model for quark phase of u, d ,s-quarks, 
because this is the first step for the study on strucures of hybrid stars with rotation and 
magnetic field in the general relativistic formulations, though it is a toy model. We use 
massless u and d quarks, and s-quark with a current mass of m s =180 MeV. We set that 
the bag constants, B, is 100 MeV fm~ 3 in this paper. Values of B > 150 MeV fm~ 3 can 
also be excluded within our model, because we do not obtain any more a phase transition 
in beta- stable matter in combination with our hadronic EOS. 

For mixed phase, we must take into account the Gibbs condition, which require the 
pressure balance and the equality of the chemical potentials between the two phases 
besides the thermal equilibrium ||21|]. We use the Thomas-Fermi-approximation for the 
density plofiles of hadrons and quarks. In each cell, we must calculate the balance of 
the colomb interaction and the surface tension. However, there are a wide range of 

uncertainties about the surface tension, o ~ 10 — 100 MeV fm~ 2 , which are suggested by 

iT — ti 

some theoretical estimates based on the MIT bag model for strangelets [22D and lattice 



guage simulations at finite temperature 11231. 12411 . 



As for this uncertainity, Maruyama et al. have showed that the surface tensions of 
a > 40 MeV fm~ 2 do not change the hardness of EOSs Q25[| . Then, we adapt two values 



of surface tension as o = 10 MeV fm~ 2 for the minimum case, and a = 40 MeV fm~ 2 
for the maximum case. 

At the maximum densities higher than two times of saturation density, muons may 
appear |2^,|23]. However, we neglect it, since the muon contribution to pressure at the 
higher density has been pointed to be very small rf28ll . 

Finally, we show our EOSs on Fig. [Q The left two panels show the pressure versus 
the baryon density for the quark-hadron mixed phase. For weak surface tension (a = 10 
MeV fm~ 2 ), the droplet structure does not appear, wheres, for strong surface ten- 
sion (a = 40 MeV fm~ 2 ), the rod structure does not appear. The number of hyperons are 
supressed because of the appearnce of quarks for both surface tensions. The right panel 
is for EOSs of nucleis with/without hyperons. 
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FIGURE 1. Baryon number density versus pressures for each matter. The left (center) panel is for an 
EOS with quark-hadron mixed phase under B= 1 00 Me V fm~ 3 and a = 1 (40) MeV fm~ 2 . The right panel 
is for hadronic EOSs with (without) hyperons shown as a solid (dashed) line. Here, we assume neutrino 
free and zero temperature for all cases; Y Ve = 0, T = 0. 

Equilibrium of magnetised rotating compact stars 

Master equations for the rotating relativistic stars containing purely toroidal magnetic 
fields are based on the assumptions summarized as follows lllql : (1) Equilibrium models 
are stationary and axisymmetric. (2) The matter source is approximated by a perfect fluid 
with infinite conductivity. (3) There is no meridional flow of the matter. (4) The equation 
of state for the matter is barotropic. (5) Magnetic axis and rotation axis are aligned. This 
barotropic condition can be maintained for our EOSs. 

NUMERICAL RESULTS 

Now, we show the configurations of strongly magnetized compact stars with/without 
quark-hadron mixed phase. In this paper, we consider the no-rotating static config- 
urations for the following two reasons. (1) Since the magnetars and the high field 
neutron stars observed so far are all slow rotators, the static models could well ap- 
proximate to such stars. (2) In the static models, one can see purely magnetic effects 
on the equilibrium properties because there is no centrifugal force and all the stel- 
lar deformation is attributed to the magnetic stress. In Fig. |2] and Fig. [3[ we show 
distributions of the baryon density and the magnetic field in the meridional planes 
for the two static equilibrium stars characterized by (1) pQ C = 1.58 x 10 15 [g/cm 3 ], 
B max = 6.2 x 10 17 G,M= 1.31Af ,/? cir = 9.97km, H/\W\ =2.11 x 1(T 3 with the quark- 
hadron phase transition (a = 40MeV/fm 2 ) [Fig. El and by (2) p 0)C = 2.56x 10 15 [g/cm 3 ], 
B max = 7.1 x 10 17 G, M = 1.31M©, # cir = 9.42km, H/\W\ = 1.64 x 1(T 3 of nucleon 
EOS with hyperons [Fig. [31 • Each physical value is as follows; the baryon mass density 
po,c, the gravitational mass M, the circumferential radius at the equator i? C7 >, the maxi- 
mum strength of the magnetic fields B max , the ratio of the magnetic energy to the gravita- 
tional energy H/\W\. These two models have a same magnetic flux of 5.00 x 10 29 G cm 2 
and a same baryon mass of 1.45M . 

Clearly, the distributions of magnetic field are different. The toroidal magnetic field 
lines behave like a rubber belt that is wrapped around the waist of the stars with hyperon 
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FIGURE 2. Distribution of (a): rest mass density [g/cm 3 ] and (b): magnetic field [G] with the quark- 
hadron phase transition((7 = 40 MeV/fm~ 2 ). 
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FIGURE 3. Same as FIG.|2]but without the phase transition. 



EOS (see Fig. [3]). However, for a hybrid star, the distribution of magnetic field has 
discontinuity for the equatorial direction (see Fig. [2]). We can understand this easily; The 
magnetic field is frozen in the matter, hence the distribution of magnetic field depend 
on the distribution of the density. Hybrid stars have discontinuities of density profiles 
because of the phase transition; e.g. in FIG.[2l the 6 km from the core of a hybrid star is 
the mixed phase matter, and the baryon density in this density region has discontinuity, 
which raise the discontinuity of magnetic field as shown in the right panel of Fig. 2. 



SUMMARY AND DISCUSSION 

In this study, we have investigated the effects of quark-hadron mixed phase on the 
magnetized rotating stars of the general relativistic equilibrium configuration. As a 
result, we find that the distribution of magnetic field for hybrid stars has a discontinuity 
in the quark-hadron mixed phase. 

Let us move on effects of strong magnetic field on EOSs. It was pointed out that 
the strong magnetic field may change EOSs Ii29ll . In their study, they estimate the 
critical strength of magnetic field which may change EOSs, ~ 10 18 G. In our study, 



we calculate the structures of magnetized compact stars around 10 18 G. Hence, such 
magnetic field may change EOSs. Especially for quark matter, it have founded that 
the energy gaps of magnetic Color-Flavor-Locked phase are oscillating functions of the 
magnetic field ||30L l31|1 . Their effects are also important, and may change our results. 

As for the temperature and lepton fraction, we assume zero-temperature and neutron 
free. We will check our result under finite temperature and neutrino trapped cases. It 
will be useful for a study on structures of proto-compact stars, or may effect on cooling 
process and cooling curves of compact stars. 
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